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Abstract. We investigate the metastability associated with the first order 
transition from normal to superfluid phases in the phase diagram of two- 
component polarised Fermi gases. We begin by detailing the dominant decay 
processes of single quasiparticlcs. Having determined the momentum thresholds 
of each process and calculated their rates, we apply this understanding to a Fermi 
sea of polarons by linking its metastability to the stability of individual polarons, 
and predicting a region of metastability for the normal partially polarised phase. 
In the limit of a single impurity, this region extends from the interaction strength 
at which a polarised phase of molecules becomes the groundstate, to the one at 
which the single quasiparticle groundstate changes character from polaronic to 
molecular. Our argument in terms of a Fermi sea of polarons naturally suggests 
their use as an experimental probe. We propose experiments to observe the 
threshold of the predicted region of metastability, the interaction strength at 
which the quasiparticle groundstate changes character, and the decay rate of 
polarons. 
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1. Introduction 

The experimental realisation of spin polarised ultracold Fermi gases has initiated 
a variety of new physics [TJ [^1 [31 HI El El El [HI El- Of particular interest is the 
understanding of strongly interacting two-component Fermi gases at zero temperature. 
Two theoretical approaches have been used to shine light on this intriguing problem: 
the study of a single atom immersed in an ideal Fermi gas of atoms in a different spin 
state and, a Monte Carlo calculation at finite polarization. The Monte Carlo approach 
revealed the theoretical phase diagram of a homogeneous Fermi g function of 

polarisation P = (JV-f — Ni)/(N^ + AQ.) and interaction strength [TU], mapping out 
a phase separation of superfluid and normal phases. The single impurity approach 
examines the quasiparticles used as the building blocks to describe these phases. In 
the strongly imbalanced limit, a single 4- fermion immersed in a Fermi sea of f fermions, 
the spin impurity atom becomes either a fermionic (polaron) or bosonic (molecule) 
quasiparticle. Complementary wave functions for each quasiparticle ([TT] [T2] and 
[131 Q31 HI] , respectively) provide groundstate energies and effective masses. Previous 
studies showed that the critical interaction strength at which the groundstate of a 
single impurity at zero momentum switches from the fermionic to bosonic branch 
[16] occurs at l/(kp^a) c ^ 0.88, with fcp-j- = (jSi^nf) 1 ^ 3 the Fermi momentum of a 
non-interacting Fermi sea of f atoms with density n-j-, and a the scattering length 
parametrizing the interaction strength between f and 4- atoms. This value is higher 
than the interaction strength ^ 0.73 at which a superfluid phase emerges in 

the limit of full polarisation |10j . An important conclusion made from this is that no 
measurement made in the groundstate will allow us to see the point at l/(kp^a) c . As 
we will show later the use of metastable and out of equilibrium processes overcomes 
this problem. 

Firstly, we determine the momentum thresholds for the decay of single impurities 
with finite momentum. We then consider the metastability associated with the normal 
to superfluid first order phase transition in the thermodynamic phase diagram. Indeed 
by applying the understanding of single quasiparticle decay to a Fermi sea of polarons, 
we predict that there exists a region where such a phase is metastable. We propose the 
use of a Fermi sea of polarons as an experimental probe to determine the threshold of 
the region. As this threshold goes to zero at l/(fc^a) c , we can observe this point for 
the first time. Beyond this threshold the decay of polarons into molecules may lead to 
a mixture of molecules and polarons, which suggests a possible way of measuring the 
molccule-polaron scattering length. In this way, we hope to open up a new regime of 
metastable physics in Fermi gases for experimental exploration. Finally, we calculate 
the decay rates for each process, within the key regions of interaction strength and 
momenta, to determine the fate of the quasiparticles. The presence of a Fermi sea of 
polarons would again be instrumental to measure the various decay rates. 

2. Background 

Our starting point is an understanding of the single quasiparticle groundstate as 
a function of going from unitarity (1/kp^a — > 0) to the "BEC" limit 

{1/k.F^a — > +oo) (see figure [l]). At l/(kp-\-a) c ~ 0.88, the critical point, the zero 
momentum energies of the polaron and molecule cross. For smaller values of 
the polaron is the groundstate of the system, and for larger values the molecule is 
the groundstate. In figure [T] and throughout this paper we use the polaron energy 
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Figure 1. Molecule and polaron zero momentum energies as a function of 
interaction strength on the "BEC" side. Red line: Molecule. Blue dashed 
line: Polaron. The dot-dashed line marks the critical interaction strength 
1/ (kpfa) c ^ 0.88. To the left of it the polaron energy is lower, and to the 
right it is the molecule that becomes the groundstate. 



calculated using the wave function proposed by Chevy [IT]. The energy in the BCS 
limit tends to Ep i = 4iran^/m, the mean-field interaction energy, and in the BEC 
limit tends to Ep i = — ^f^- Here and in the following, we take H = 1. The 

effective mass of a polaron (m Pol ) obtained from |15j becomes just the bare mass 
in the BCS limit and diverges at l/kp^a ~ 1.17. For the molecule, energies and 
effective masses are obtained from jT3l [14j [15] . In the BEC limit the energy tends to 
EmoI = — £Ff' The effective mass equals the bare molecule mass ra* Mol = 2m 

in that limit, while approaching unitarity it grows, and it diverges at 1/kp^a ~ 0.55. 
These masses and energies for both polarons and molecules have been found to be 
very accurate by comparison with Quantum Monte Carlo calculations [101 116] . 

3. Decay processes and thresholds 

We begin by considering the decays that the polarons and molecules can undergo 
when they are not in the groundstate (e.g. at finite momentum). The decay process 
of a quasiparticle can only occur when it satisfies energy and momentum conservation. 
This is generally only possible when the initial quasiparticle momentum p is above a 
threshold momentum Pph- 



Metastability in spin polarised Fermi gases and quasiparticle decays 



4 



3.1. Polaron decay 

For the polaron at zero temperature, the decay processes we consider are: 

A. Polaron — > Molecule + hole 

B. Polaron — > Molecule + 2 holes + particle 

A finite momentum polaron also has a finite relaxation time as described in |17| , which 
schematially reads 

C. Polaron (p) — > Polaron {p' < p) + hole + particle 

The three-body decay of a polaron into a molecule, two holes and a particle 
has been considered in [18] in the special case p = 0. Generalizing this calculation 
to non-zero momentum, we find that process B can occur at any momentum when 
AE = Ep i — Em i > 0, as expected, and only for 

p > P* h = y/-2m* Pol AE (1) 

when AE < 0. Here Ep Q i and EmoI are the energies of a zero-momentum polaron 
and molecule at the given value of interaction strength. 

For process A, the two-body decay of a polaron with momentum p into a molecule 
and a hole, conservation of energy and momentum require the following equality to 
be verified: 

p 2 (p + q) 2 

Epot + 2^ + ^ = Emo[ + 2M* ' (2) 

2 

where £ g = — ep^ is the kinetic energy of a majority particle measured with respect 
to the Fermi surface. The minimum momentum Pp h at which process A is allowed 
is found by setting the hole at the Fermi surface, and by taking its momentum to be 
anti-parallel to the one of the polaron: 

P Th 2 , rp _ ( P Th ~Pf) , , 

—+E Po i- — — +E M oi- (3) 

The processes including more particle-hole pairs (e.g. Molecule + 3 holes + 
2 particles) have the same energy threshold as the state resulting from process B, but 
lower rates so we do not consider them here. 

In conclusion, polarons with non-zero momentum are stable towards A/B decay 
for momenta smaller than P^ B . 

Figure [2] shows our results for the polaron decay thresholds, as given by ([T]) 
and In the region below l/(fc^a) c , where the zero momentum polaron is the 
groundstate of the system, a polaron at finite momentum can nevertheless be unstable 
to decay processes A, B, and C. The solid blue line gives the momentum threshold 
P^ h for the polaron to decay via process B. On this line, decay results in a zero 
momentum molecule for 1/ (fc^a < \/{kp^a) c . Finite momentum molecules also result 
from process B above the solid blue line. The blue dashed line gives the momentum 
threshold Pjt h for process A, which generally creates a molecule at finite momentum 
even on the threshold. Decay processes including more particle-hole pairs all set in 
at momenta above the solid blue line, i.e. for p > P^ h . For l/[kp^a > l/(kp^-a) c , a 
polaron at zero momentum is unstable to a process B decay into a molecule at finite 
momentum whereas process A continues to affect only higher momentum polarons. 
Polarons with any finite momentum are unstable to momentum relaxation (process C) 
on both sides of l/(fcj?fa) c . The momentum thresholds for processes A and B become 
equal where the molecule's effective mass diverges (1/kp^a ~ 0.55). 
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3.2. Molecule decay 

The stability of a molecule is calculated in a similar way. The decay channels we 
consider are: 

A. Molecule Polaron + particle 

B. Molecule —¥ Polaron + 2 particles + hole 

C. Molecule (p) — > Molecule (p' < p) + hole + particle 

For 1/kpa < l/(/cFt a )c, a, zero momentum molecule decays via process B to a polaron 
with finite momentum, and via process A at higher momentum. For 1/kpa > l/kp^a c , 
process B preceeds process A with increasing momentum until where the polaron's 
effective mass diverges. In analogy with the polaron decay considered above, the 
momcnum thresholds for the molecule are determined by energy and momentum 
conservation. They are, 

P Th ( P Th-Pp) 
'■" 1 E Mo i = ?r— \-Epoi, (4) 



1m* 1m 



^^(M? (5) 

where P^ h — for l/kp^a ^ l/(kp^a) c . Again, the threshold momentum for process 
A is found when the majority particle is on the Fermi surface. 

Our results for the threshold momenta for polaron and molecule decay, as given 
by 0, Q and, (JBJ are shown in figure [2j 



4. Metastability of polaron gas 

Thus far we have analysed the decay processes of single quasiparticles. The single 
quasiparticle problem is a limiting case of the spin imbalanced Fermi gas. The 
groundstate phase diagram of a spin imbalanced Fermi gas, with polarisation versus 
interaction strength, has been calculated in [1IT| using a Monte Carlo approach. In 
this calculation, the J, atom concentration is kept finite for a macroscopic system even 
in the P — > 1 limit. We now examine how the stability analysis above is connected to 
this equilibrium phase diagram. 

Our particular consideration is the normal to superfluid first order phase 
transition predicted by the Monte Carlo calculations [TO]. This transition from a 
partially polarised normal phase to a state with separated superfluid and normal 
phases is represented by the dashed green line in figure [3| To make this line in figure 
|3j we have assumed that the polarons form a weakly interacting Fermi sea, so that we 
can convert the critical density at which phase separation occurs, into a critical Fermi 
momentum pp° l (l/kp^a) for the polarons. In the phase separated state, the Monte 
Carlo calculation includes the interactions between the molecules in the condensate. 
In contrast, in the single J, atom calculation there is at most one molecule and therefore 
no condensate. 

If we assume that the relevant processes in which the polarons can be converted 
into molecules, within experimentally realistic timescales, are only the ones considered 
in section |3.1| and that these single impurity processes can be used to analyse the 
stability of the Fermi sea of polarons (neglecting for instance multiple polaron decay) , 
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Figure 2. Momentum thresholds for various decay processes for both polarons 
and molecules about the critical point. Blue solid line: Polaron — > Molecule + 
2 holes + particle; Blue dashed line: Polaron — > Molecule + hole; Red solid 
line: Molecule —> Polaron + 2 particles + hole; Red dashed line: Molecule — > 
Polaron + particle. The groundstate of the system changes from being a polaron 
to a molecule at l/kp^a ~ 0.88, beyond which even a polaron at p = can 
decay via process B; otherwise, only polarons at p > are unstable in this way. 
Process A is relevant at higher polaron momentum and generally results in a finite 
momentum molecule even on the threshold. The process A and B thresholds meet 
when the molecule's (polaron's) effective mass diverges at 1/kp^a ~ 0.55(1.17). 
Note that for polarons at large p (~ PFf) > momentum relaxation processes are 
very strong so that such a polaron is no longer a well-defined quasiparticle. 



this gives rise to a region of metastability in the phase diagram. If a Fermi sea of 
polarons prepared in the groundstate (below the green dashed line) is adiabatically 
taken above the dashed green line by increasing 1/kpa, we expect it to persist as 
a metastable state. Even though the density of the \, atoms is so high that the 
true groundstate is a phase separated state, the state is stable since there are no 
polarons with large enough momentum to decay to a molecule via process B. The 
decay of polarons to molecules sets in only when l/kpa is increased further so that 
the Fermi momentum of the polarons is larger than the momentum threshold for 
process B (pp° l = Pp^)', we ig nore here molecule-polaron interactions (see below). 
This momentum is given by the solid blue line in figure [3j At this point, the highest 
momentum polarons will decay into molecules. 

This motivates the following proposal to use a Fermi sea of polarons as an 
experimental probe to explore the metastable region and observe molecule formation 
through process B. Due to the Pauli exclusion principle, the momentum relaxation 



Metastability in spin polarised Fermi gases and quasiparticle decays 



7 



process (C) is suppressed allowing for an analysis of the other decay processes. Let 
us suppose for the moment that the molecules into which polarons can decay do not 
interact with the polarons. Then, to determine the shape of Pp h , we can prepare 
a polaron Fermi sea with a Fermi momentum smaller than the dashed green line 
in figure [3] at unitarity, for example. Then we increase adiabatically. The 

system will cross the first order phase transition line, but since it is metastable (as 
discussed above), it will remain a Fermi sea of polarons rather than form a condensate 
of molecules. As we continue beyond the phase transition, the Fermi momentum will 
become equal to the momentum threshold for process B (pp° l = Pph)- Increasing 
infinitesimally beyond this point will now lead to the decay of polarons at 
Pp° l . The polarons on the Fermi surface decay into molecules at zero momentum. 
The appearance of these molecules can then be detected experimentally as the tell-tale 
sign of the solid blue line Pp h - The experiment can be repeated using an appropriate 
Pp° l to find P® h at different interaction strengths, making sure that only a small 
number of molecules are created each time (so as to be able to ignore effects beyond 
the threshold, e.g. molecule-molecule interactions) but a large enough number to be 
observable. The size of the polaron Fermi sea needed will decrease as P? h decreases, 
so that the value l/(kp^a) c is found in the limit of a single polaron. 

If we now take into account molecule-polaron interactions the value of pp° l 
at which molecules are first formed will change. If we consider one molecule in 
the final state, the (unknown) molecule-polaron interaction changes the energy by 
AE — gpMnpoi (assuming a mean-field approximaton). This will lead to a positive or 
negative shift of the threshold curve of production of molecules (see dot-dashed blue 
lines in figure [3]). This shift tends to zero at l/(k F ^a) c where n Poi —> and where 
the threshold curve meets Pp h - One could in principle use the difference between the 
experimentally observed curve and P^ h (which is known theoretically) to determine 
this shift and so the molecule-polaron scattering length. Note that we have ignored 
the effects of the polaron-polaron interaction, which are known to be small. 

We have therefore established that the region between the dashed green and dot- 
dashed blue lines in figure [3] represents a metastable phase consisting of a Fermi sea of 
polarons. Such a metastable phase is characteristic of a first order phase transition. 
It is important to note that the state may well be metastable beyond the dot-dashed 
blue line since an analysis of its stability in that region would require us to take into 
account the presence of a finite quantity of molecules. We also raise the intriguing 
possibility of a final state containing the remaining Fermi sea of polarons coexisting 
with a condensate of molecules, within a background of f particles. 

5. Single quasiparticle decay rates and experimental observability 

5.1. Decay rates calculation 

The phase diagram is animated by considering the rates of each decay process. The 
decay rates of zero momentum polarons and molecules via process B are presented 
in |18j . The decay rate of a polaron with momentum p through process A is given 
by the imaginary part of the on-shell polaron self-energy £p oJ shown in figure |4^i, i.e. 

TpoZ (p) = — ImSp j (p, Ep i + — ^ ■ At zero temperature, one obtains 

rffi (p) = - W / s(ae + + e, - , (6) 

Jq<k F ^ \ Zm Pol Zm Mol ) 
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Figure 3. Stability of Fermi sea of polarons: The dashed green line shows the 
first order transition from normal to phases including a superfluid (Fig. 4 in 
I10|). The solid blue line is the Fermi momentum of a Fermi sea of polarons equal 
to the momentum threshold for process B (pp° l = P^ h )- The dot-dashed blue 
lines represent how the solid blue line might be shifted due to molecule-polaron 
interactions. 



which is a simple Golden rule expression. Here, g = — y is the atom-molecule 

coupling in vacuum |18j . and m r is the atom-molecule reduced mass. To derive this 
equation, we have used a pole expansion of the molecule propagator with quasiparticle 
residue ZmoI ■ The decay rate may be calculated analytically. Assuming for simplicity 
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Mol 



2m-f and rn* Pol 



m-f, we obtain 
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Tpoz (p) = Z M ol £Ft 



IP- 



pA \ 
r Th) 



(Pf 



Th 



+ 2k Ft - p) 



for P£ h <p < Pfi h 



process is PA = k 



2kp-f and zero otherwise 
1 



pk Ft ' ^ 

The momentum threshold for this 



2 1- 



AE 



Note that Pfi h > even when AE > 0, 



and that one has < 0.5 for all scattering lengths [T51 116| . 

The decay rate of a molecule with momentum p via process A, the creation of 
a polaron and a majority particle, _i_s given by the imaginary part of the molecule 

A calculation analogous to the polaron case 



self-energy S^ /oZ depicted in figure |4 
considered above yields, 



p-pA 



p-pA 



4k 



Ft 



kp^a 



pk F ^ 



(8) 
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(b) Molecule self-energy 



y(B) = 

^Mol 





Figure 4. Diagrams showing the decay processes A an B for (a) polarons and (b) 
molecules. The thin lines represent majority atoms, wavy lines a molecule, and 
thick lines a polaron. The polaron-molecule matrix element g is represented by 
a thick dot and the square represents off-resonant scattering between a polaron 
and the majority atoms. 



with Zp i the quasiparticle residue for the polaron. The threshold momentum for this 

fi 1 



process is P, 



A 
Th 



t F f J 



with T Mol 



for p < P 



A 

Th- 



Note that 



in this case, as opposed to the polaron case, there is no maximum momentum for 
the molecule above which there is no decay via process A. This is because for large 
momentum, the molecule can always dispose its energy and momentum to a majority 
particle, whereas the polaron has to dispose it into a majority hole within the Fermi 
sea. 

The corresponding self-energy for process B, Sp o; is shown in figure [4^,. 
Generalising the calculations in [18] to non-zero momentum p gives 



r Poi (P) °c Z Mo ieFt 



[AE 



2m; 



7.) 



(9) 



The analysis in [18] showed that the proportionality constant in this expression is of 
order unity. 

Likewise, the decay rate of a molecule via process B can be obtained by 
generalising the calculations in |18j to a finite momentum p. We obtain 



T Moi (?) « Z Po ie F1 - 



-AE- 



p_ 

2m , 



(10) 
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where the constant of proportionality is, again, of order unity. P^ h {P^h) ' s the point 
at which Tp ol (Tf /oi ) equals zero. 

For a finite momentum, the polaron can scatter off the majority particles giving 
rise to momentum relaxation with a rate — ^— (process C). The high velocity regime 
kFi "C nijV -C kp\ analysed in [XT] determines the rate of relaxation of a single 
impurity. If we consider only a single impurity, spin statistics becomes redundant, 
allowing the same calculation to be used for the molecule. The momentum relaxation 
rate, at T = for molecule or polaron is then, 

1 2,7T 2 m Mol/Pol vA , . 

= ^n\ — jk ' ( n ) 

T Mol/Pol «JO K Ft 

where 7 is determined by the scattering amplitude U — 8< *'^ /F ° l = — 7. At 
unitarity, we have [ip i = —aep^ with a ~ 0.6 [15]. One can also use HmoI = 
1 + 5™^ valid for -r^— > 0.7 Ql] (where a ~ 1.18a) to find the molecule 



momentum relaxation rate in terms of the interaction strength, 

— = ^(^) 2m f^. (12) 
TMol "On kf^ 

The high power of the velocity v and effective mass m* indicate impurities at large 
momenta p ~ pp-\ are no longer well-defined quasiparticles. 

In figure [5] we plot the decay rate of polarons and molecules via processes A and B 
as a function of momentum p for various AE. We see that once p > P^ h and process A 
sets in, it quickly dominates over B. In fact, once active, process A dominates over B by 

several orders of magnitude since the ratio is, to lowest order oc ( p ~ p th \ (2^t\ , 

6 ' r « V Pft / \PThJ 

This is expected since process A is a 2-body process, whereas process B is a 3-body 
process. On the other hand, for p < P^ h process B of course dominates as process 
A is not allowed. However, since the typical time scale for process B is very long, of 
order 10-100ms, the momentum relaxation of the polaron via process C is in general 
much faster. This justifies our proposed experiment to determine P^ h using a Fermi 
sea of polarons in equilibirum, to suppress process C. 



5.2. Decay rates experiment 

The polaron Fermi sea can also be used to measure the rate of process B. The initial 
state is a polaron Fermi sea prepared in the groundstate (i.e. below the dashed green 
line in figure [3]). We then instantaneously increase the interaction strength so that the 
some fraction of the polarons have a momentum p > P^ h and so are then susceptible 
to decay via process B. Alternatively, by increasing the interaction strength to above 
l/(kFfa) c , every polaron has a momentum p > P^ h , and so the whole Fermi sea 
is susceptible to decay via process B. The molecules resulting from this decay are 
unable to decay via process A or £?, however they lose momentum via process C . The 
resultant state is therefore expected to be a condensate of molecules. A measurement 
of the initial growth rate of the number of molecules or the loss rate of polarons 
determines the rate of process B for polarons, averaged over the momenta of the 
polarons decaying, at a given l/fc^a. 

Similarly, the rate of process A can be measured, with some fraction of the 
polarons having a momentum p > P^ h . In this case, the Fermi sea of polarons 
will decay via both processes A and B. The two processes can be distinguished since 
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Figure 5. Decay rates of finite momentum polarons (top) and molecules (bottom) 
for processes A (continuous line) and B (dashed line). From left to right, we have 
AE/ep^ = —0.2, 0, 0.2. The thick dot indicates the threshold momentum for 
process B, P^ h - Process A, above its threshold dominates over B for p < PFf- 



process A has a significantly faster rate and typically results in molecules at finite 
momentum. 

6. Conclusions 

In this paper, we have shown how quasiparticles in a two-component Fermi gas behave 
at finite momenta in the limit of extreme imbalance. By considering energy and 
momentum conservation of single impurities, we determined the momentum thresholds 
beyond which the quasiparticles are susceptible to the most significant decay processes, 
and we calculated the rates of each. Using this and assuming we can use single 
quasiparticle decay processes to describe the decay of a Fermi sea of polaron, we have 
identified a region of metastability for the partially polarised normal phase about the 
normal to supcrfluid first order transition. We then described how the Fermi sea can 
be used as an experimental probe to observe the metastable region, determine the 
crossing point of single impurity energies 1/ (kF^a) c for the first time and measure the 
decay rate of process A and B. The phase diagram we construct and the adiabatic 
experiment we propose to explore it can be used as a measure of the unknown molecule- 
polaron interaction strength, and leaves open the possibility of observing a novel state; 
a mixture of molecules and polarons. 
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